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The  deaign  of  correlation  regions  for  tracking  with  cooperative  uaera  ia  examined 
under  the  requireoMnt  to  design  the  correlation  regiona  on  a conatant  probability 
of  intercept  baaia.  Starting  vith  the  aaaumption  of  independent  and  Gausaian- 
diatributed  range  and\aaiakuth  errora  in  the  aenaor  and  assuming  a conatant- 
coefficient  iaotropic  d -S  tracking  filter,  it  ia  ahovn  hov  the  correlation  region 
deaign  muat  include  auch  factora  aa  aenaor  errora,  timing  jitter,  tracking  errora, 
and  the  aaynchronoua  operation  of  the  tracking  function  vith  reapect  to  the  aen- 
aor meaaurementa . Uaing  a apecific  example  to  obtain  numerical  reaulta,  it  ia 
ahovn  that,  in  the  caae  of  atraight-line  trajectoriea , the  radiua  of  a circular 
correlation  region  ia  linearly  proportional  to  the  distance  from  the  aenaor  in 
the  region  vhere  aximuthal  errora  predominate  and  ia  conatant  in  the  region  vhere 
range  errora  predominate.  For  a maneuvering  target,  it  ia  ahovn  that  the  aice  of 
the  correlation  region  muat  be  equal  to  the  aum  of  the  radiua  used  for  the 
atraight-line  caae  plua  the  magnitude  of  any  tracking  biaa  which  reaulta  because 
of  the  deviation  from  a atraight-line  trajectory  aa  aaaumed  in  the  tracking 
filter.  By  examining  various  types  of  maneuvers,  an  upper  bound  Is  derived  for  the 
magnitude  of  the  biaa  vhich  could  reaaonably  be  expected  in  typical  maneuvera. 

By  apecifying  the  aise  of  the  correlation  region  on  a conatant  probability  biaa, 
it  ahould  be  poaaible  to  obtain  better  diacrimination  against  false  targets  and 
improved  detection  of  maneuvers  by  sensing  the  development  of  tracking  biases. 
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EXECUTIVE  SUMMARY 


The  baeie  for  all  advanced  air  traffic  control  functions  ia  the  ability  to 
predict  the  position  of  an  aircraft  sufficiently  far  into  the  future  to  allow 
intervention  by  a controller  in  situations  in  which  this  is  warranted.  The 
prediction  of  future  position  is  based  on  velocity  estisMtes  obtained  from 
a tracking  filter  which  estisMtes,  via  nuawrical  differentiation,  the  time 
derivatives  of  the  position  reports  for  a given  aircraft.  Within  the  tracking 
algorithm,  the  correlation  region,  or  search  area,  is  one  of  the  means  by  which 
new  target  reports  are  associated  with  exiiting  tracks  derived  from  previous 
reports.  The  topic  of  interest  in  this  paper  is  the  design  of  search  areas 
which  will  achiave  a high  level  of  correct  association  between  the  existing 
tracks  and  new  target  reports,  but  which  are  sufficiently  selective  to  reject 
erroneous  reports  such  as  false  targets.  In  order  to  achieve  this  goal  it 
ia  necessary  to  design  search  areas  which  are  adaptive  to  the  statistical 
performance  throughout  the  coverage  area  of  the  sensor. 

The  design  of  variable  or  dynamic  search  areas  ia  based  on  both  the  measure- 
ment errors  of  the  sensor  and  tha  performance  of  the  tracking  filter  for 
various  types  of  maneuvers.  In  all  practical  tracking  algorithms  the  trajec- 
tories of  maneuvering  targets  are  more  difficult  to  follow  than  straight- 
line  tracks.  Since  the  sise  of  the  dynamic  search  area  varies  with  the 
sice  of  the  expected  measurement  error,  and  other  factors,  this  technique  is 
able  to  provide  a batter  discrimination  between  straight-line  and  maneuvering 
targets  which  results  in  an  improvement  of  the  tracking  performance  (via 
a reduction  in  the  bias  observed  in  heading  for  maneuvers  and  a reduction 
in  the  heading  jitter  for  straight-line  tracks).  Since  the  search  area 
design  is  based  on  known  statistical  characteristics  of  the  sensor  measuresient 
errors,  the  theoretical  performance  of  the  tracking  algorithm  and  the  devi- 
ations from  the  theoretical  performance  encountered  in  practical  tracking 
algorithms,  the  dynamic  search  area  concept  guarantees  that  a specified 
level  of  tracking  perforMnce  will  be  maintained  throughout  the  coverage  area 
of  the  sensor.  The  numerical  techniques  developed  in  this  paper  are  appli- 
cable to  the  design  of  dynamic  search  areas  for  use  with  tracking  filters 
using  either  DABS  or  ATCRBS  data.  In  a subsequent  paper,  the  specific  levels 
of  tracking  perforMnce  improvement  which  can  be  obtained  using  the  dynamic 
search  area  will  be  found. 
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1. 


INTRODUCTION 


The  process  of  radar  data  correlation,  which  is  defined  as  the  association 
of  surveillance  data  with  a track  (a  mathematical  estimation  of  a target 
trajectory),  has  been  considered  previously  [1-7],  The  requirement  for 
a correlation  process  usually  arises  in  the  transfer  of  tracking  respon- 
sibility between  surveillance  volumes  covered  by  different  sensors  or  in 
the  multiple-target/multiple-sensor  case  for  the  same  surveillance  volume 
where  an  association  between  tracking  and  surveillance  data,  which  may  contain 
false  reports,  must  be  made.  In  most  cases  the  size  of  the  gate  (or  search 
area)  in  which  correlation  will  be  allowed  is  based  on  the  statistical  prop- 
erties of  the  residual  errors  between  the  measured  target  data  and  the  pre- 
dicted state  of  the  trajectory  at  the  time  of  interest.  If  it  is  assumed 
that  the  distribution  of  the  residuals  is  Gaussian  and  independent  in  each 
dimension,  then  the  gate  boundary  will  usually  be  taken  as  an  equiprobability 
contour  with  the  statistical  properties  based  on  the  chi-square  distribution 
(3,  5,  and  7]. 

The  usual  context  for  correlation  analyses  is  in  the  case  of  uncooperative 
users  in  which  there  is  a noise  background  of  spurious  data  combined  with 
true  target  reports  corrupted  by  measurement  errors.  For  the  purpose  of 
this  study,  however,  it  will  be  assumed  that  cooperative  targets  (i.e., 
transponder  equipped)  are  being  observed  and  for  this  the  correlation  prob- 
lem is  of  a slightly  different  nature.  The  system  of  interest  in  this  case 
is  the  Air  Traffic  Control  Radar  Beacon  System  (ATCRBS)  [8]  in  which  the 
correlation  process  is  used  to  reject  data  which  are  unreasonable  because 
of  inherent  system  limitations  (e.g.,  identification  errors  in  the  beacon  code 
reply  can  be  caused  by  overlapping  replies  from  closely  spaced  targets,  erro- 
neous positions  can  be  reported  due  to  reflections,  etc.).  In  addition  to 
rejecting  unreasonable  data,  the  correlation  process  must  also  determine  which 
of  several  alternative  data  points,  if  any  exist,  is  the  best  for  a particular 
track  and,  also,  to  detect  a target  maneuver  by  the  size  of  the  residual 
error.  It  has  been  suggested  that  data  association  and  tracking  should  be 
treated  in  a unified  manner  rather  than  relying  on  simulation  of  alternative 
techniques  [1],  and  it  is  the  objective  in  this  study  to  illustrate  the 
interrelationship  between  correlation  and  tracking  and  to  show  how  certain 
practical  problems  of  implementation  can  affect  the  design. 


2. 


MATHEMATICAL  DEVELOPMENT 


2.1  SURVEILLANCE  MODEL. 

In  order  to  make  probabilistic  statements  about  the  size  of  the  correlation 
regions,  some  form  of  statistical  model  must  be  used  to  describe  the  sensor 
measurement  process.  It  will  be  assumed,  as  it  has  been  in  other  studies 
[20],  that  the  surveillance  system  measures  target  position  in  a polar 
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coordinate  system  and  that  the  range  and  atimuth  errors  are  small  in  magni- 
tude, Gaussian  distributed,  and  statistically  independent.  The  transformation 
to  the  Cartesian  coordinate  system  used  for  tracking  is,  where  8 is  measured 
from  the  y-axis, 


x * p sin  8 

(1) 

y • 0 cos  8 

(2) 

and  for  small  errors  a first-order  approxiaiat ion  is 

Ax  • sin  8 A p ♦ p cos  8 A0 

(3) 

Ay"  cos  8 Ap-psin0A0. 

(4) 

It  is  implicitly  assumed  in  the  above  transformation  that  data  from  targets 
at  different  altitudes  have  been  mapped  onto  a common  plane  using  an  appro- 
priate projection  technique.  Assuming  that  the  errors  in  range  and  asimuth 
are  unbiased  then  Ax  and  Ay  are  bivariate  Gaussian  random  variables  with 

xero  swan  and  variances 

02  “ sin  20  ♦ ( P cos  0>2 

(5) 

Oy  • cos  2 8 Of?  ♦ (p  sin0)2  Oq 

(6) 

2 2 2 

oXy  • (ap  - P 0 q)  sin  8 cos  8 • 

(7) 

Equations  (3)  and  (4)  can  be  solved  simultaneously  for  A P and  A $,  and  the 
resulting  equations  can  be  used  to  express  the  variances  in  a polar  coor- 
dinate system  in  terms  of  the  variances  in  a Cartesian  coordinate  system; 
i .e. , 


a 

P 

a 2 
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0 2 sin20  ♦ a 2 cos20  ♦ 2o 
x y xy 


sin  $ cos  8 


(o' 


2 2 
cos  8 * ° 

y 


sin20 


2 o sin  8 cos  8 )/p  • 
xy 


(8) 


(9) 


While  (8)  and  (9)  are  valid  for  any  Gaussian  distribution,  it  will  not  neces- 
sarily be  true  that  P and  8 are  statistically  independent.  If,  however,  the 
statistical  characteristics  in  the  Cartesian  coordinate  system  were  originally 
derived  from  a polar  coordinate  system,  via  (5)  to  (7),  then  p and  8,  in  a 
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■ polar  coordinate  system  with  statistical  characteristics  defined  by  (8)  and 
(9),  will  be  statistically  independent.  The  usefulness  of  the  t rans format  ion 
from  a polar  to  a Cartesian  coordinate  system  and  subsequently  reverting  to  a 
polar  coordinate  system  will  become  apparent  in  a later  section  in  which  the 
evaluation  of  certain  integrals  is  discussed. 

2.2  SPECIFICATION  OF  THE  TRACKING  FILTER. 

The  sensor  measurements  consist  of  the  sequences  { p(k)}  and  (0(k)}  where  p(k) 
and  0(k)  are  the  true  range  and  azimuth,  respectively,  of  the  target  at  time 
epoch  k.  In  practice,  the  sensor  measurements  are  corrupted  by  errors,  Ap(k) 
and  A0( k),  and  it  will  be  assumed  that  the  sequences  (Ao(k)}  and  (A0(k))  are 
each  white  and  stochast ical ly  independent  of  one  another.  After  t rans format  ion 
to  a Cartesian  coordinate  system  via  (1)  and  (2),  the  measurement  sequences 
(X_(k)}  and  {Yn(k))  with  the  associated  error  sequences  as  defined  by  (3)  and 
(4)  are  input  to  a digital  filter  which  is  characterized  by  the  recursive 
equat ions , 

X (k)  - X (k)  ♦ •»  (X  (k)  - Xn(k) ) (10) 

s P m ” 

V (k)  - V (k-i)  ♦ (6/T)  (X  (k)  - XD(k))  (11) 

s s m P 

Xp(k*l)  - X (k)  ♦ T V (k)  (12) 

v s s 

where  XB,  VB,  Xp  are  the  smoothed  position,  smoothed  velocity,  and  predicted 
position,  respectively,  and  similarly  for  Y.  The  sensor  data,  X^,,  are  assumed 
to  be  available  at  a constant  rate  specified  by  the  time  interval  T.  The 
properties  of  this  filter,  known  as  ano-8  filter,  are  well  known  19-24).  The 
usefulness  of  this  filter,  as  compared  to  others  with  superior  performance, 
lies  mainly  in  its  ease  of  implementation  and  limited  computational  require- 
ments which  allows  tracking  of  large  numbers  of  targets  with  computers  of 
relatively  modest  capacity.  The  smoothing  constants  <•  and  8 determine  the 
performance  of  the  tracking  filter  and  various  criteria  have  been  developed  for 
specifying  these  constants  (9-13).  Smaller  smoothing  constants  generally  give 
better  noise  reduction  properties,  while  larger  smoothing  constants  give  better 
transient  or  maneuver-following  capability.  Since  it  is  not  possible  to  select 
smoothing  constants  which  are  optimal  in  all  cases,  it  is  frequently  necessary 
to  use  several  seta  of  smoothing  constants  to  achieve  a practical  system.  The 
choice  between  sets  of  smoothing  constants  is  then  based  on  the  magnitude  of 
the  difference  between  the  measured  position  and  the  predicted  position  which 
is  referred  to  as  the  track  datum  deviation. 

In  the  analysis  of  the  >»- 8 tracking  filter,  it  is  normally  assumed  that  the 
filter  operates  at  a constant  rate  with  data  available  at  the  same  rate. 

While  this  is  a reasonable  assumption  in  the  case  of  a single  sensor,  in  the 
case  of  multiple  sensors,  each  of  which  covers  only  a portion  of  the  total 
surveillance  volume,  it  is  unlikely  that  the  sensor  system  and  the  tracking 
system  will  operate  in  a synchronous  manner  so  that  the  t ime  of  receipt  of  the 
data  may  not  be  the  same  as  the  prediction  time  of  the  tracking  system.  Such 
a situation  arises  in  the  enroute  air  traffic  control  system  in  which  the 
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•mooching  and  prediction  process  is  done  ac  fixed  intervals  which  are  asyn- 
chronous with  respect  to  the  sampling  intervals  of  the  sensors  (in  fact,  the 
sampling  intervals  of  the  sensors  are  not  all  the  same).  To  compensate  for 
the  asynchronous  operation  of  the  tracking  and  sensor  systems,  the  estimated 
velocity  is  used  to  shift  the  measured  position  in  time  so  as  to  be  time 
coincident  with  the  predicted  position,  thus  attempting  to  duplicate  the 
operation  of  a synchronous  sensor/tracking  system.  In  this  case,  (10)  and 
(11)  become 

Xa(k)  - Xp(k)  ♦ <>  (Xgfk)  ♦ ATc(k)  V,(k-1)  - Xp(k))  (13) 

V8(k)  - V,(k-1)  ♦ ( B /T)  (X,„(k)  ♦ ATc(k)  V,(k-l)  - Xp(k))  (14) 

where  ATC  is  the  difference  between  the  time  of  receipt  of  Xgfk)  and  the 
time  used  in  the  calculation  of  the  predicted  position  Xp(k).  The  require- 
ment to  adjust  the  smoothing  and  prediction  process  to  account  for  the  time  of 
actual  receipt  of  the  datum  has  been  recognized  along  with  the  fact  that  the 
gate  size  must  be,  "...large  enough  to  account  for  all  sources  of  error"  (10] 
in  order  to  achieve  correlation  with  the  target  data. 

The  fact  that  a 'time-correction'  process  must  be  used  to  compensate  for  the 
asynchronous  operation  of  the  tracking  algorithm  implies  that  the  correlation 
process  must  be  a two-stage  procedure,  since  the  process  of  time  correction 
itself  implies  that  at  least  a tentative  correlation  has  already  been  made, 
unless,  of  course,  an  exhaustive  searcn  is  made  considering  every  possible 
track-datum  pairing.  Many  techniques  have  been  developed  for  making  tentative 
correlations  and  these  include:  (a)  associating  the  beacon  code  of  the  datum 
with  the  code  of  a track  which  already  exists,  (b)  dividing  the  surveillance 
volume  into  subvolumes  with  lists  of  the  tracks  in  each  subvolume,  and  (c) 
the  use  of  a preliminary  search  area  defined,  for  example,  as 


^(k)  - Xp(k) j < D 

(15) 

|Ym(k)  - \p(k)|  < D 

(16) 

where  D is  chosen  large  enough  to  include  all  possible  track-datum  pairings 
which  might  be  allowable  in  the  second  stage  of  correlation.  It  will  be 
assumed  that  some  technique  is  available  to  make  tentative  correlations, 
similar  to  the  coarse/fine  correlation  in  [3]  or  one  of  the  techniques  above, 
so  that  it  is  not  necessary  to  make  exhaustive  search  of  every  possible  track- 
datum  pairing.  Regardless  of  which  technique  is  used  the  characteristics  are 
the  same,  namely  that  the  preliminary  or  coarse  stage  must  process  large 
amounts  of  data  with  very  low  computational  requirements  for  each  pairing, 
while  the  second  or  "fine"  stage  processes  small  amounts  of  data,  but  may  have 
considerably  higher  computational  requirements  for  each  pairing. 
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2.3  TRACK  DATUM  DEVIATION. 


In  Che  aecond  or  fine  stage  of  the  correlation  process,  which  is  of  primer) 
interest  in  this  study,  some  means  must  be  available  to  assess  the  quality 
the  various  track/datum  pairings  which  result  from  the  tentative  correlatic 
made  in  the  first  stage  of  the  correlation  process.  Assuming  that  all 
other  relevant  factors  are  equal  the  best  track/datum  pairing  will  be  that 
which  minimises  the  magnitude  of  the  track  datum  deviation,  A" r",  which  is 
given  by  (including  the  time-correction  factor), 

A"r(k)  -^(k)  ♦ ATc(k)vs(k-l)  - rp(k)  (17) 

where 'r‘m,  ^Ts,  and  Yp  are  the  measured  position,  smoothed  velocity,  and 
predicted  position,  respectively,  and  it  is  now  explicitly  denoted  that  the 
track  datum  deviation  is  a vector  as  are  the  positions  and  velocity.  Since 
|£r(k)|  will  be  used  as  the  final  basis  for  correlation,  the  statistical 
characteristics  of  &~r(k)  must  be  determined  to  know  how  large  a deviation  is 
reasonable.  Expressing  each  of  the  terms  in  (17)  as  the  sum  of  a true  com- 
ponent plus  a random  error  yields  (where,  for  simplicity,  the  dependence  on 
k has  been  assumed), 

hr  ■ "r^j  ♦ "r’ei  + (ATf+O  (Vj+Ve)  - "r^  “~re2  (18) 


or 


Ar*  * rjj  ♦ ATfVT  - r^^  + ♦ t Vf  ♦ ATjVe  ♦ c Ve  - 


(19) 


where  < is  a timing  error  which  arises  from  the  fact  that  time  is  measured  in 
discrete  increments.  Of  course,  in  a digital  computer  all  data  are  quantized 
but  the  additional  errors  introduced  by  finite  precision  computations  will 
be  assumed  to  be  negligible.  If  the  trajectory  dynamics  implied  by  (12), 
i.e.,  a constant  velocity  straight-line  path,  are,  in  fact,  correct,  then  the 
first  three  terms  of  (19)  will  vanish  since  Tj+  Tt“t2-  If  the  target 
is  maneuvering,  however,  there  will  be  a bias  error  introduced  in  Ar  by  the 
fact  that  the  assumed  system  model  is  not  correct.  For  the  present,  however, 
it  will  be  assumed  that  the  system  model  is  correct  so  that  (19)  becomes 

A"?  ■ "r^j  ♦ f Vf  ♦ TfVe  * *^e  “ re2  (20) 


in  which  the  statistical  properties  of  r^  are  given  by  (5)-(7).  In  the 
case  of  the  timing  error,  ( , it  will  be  assumed  that  these  errors  are  uni- 
formly distributed  over  an  interval  ATq  with  E(O*0,  which  can  easily  be 
obtained  by  an  appropriate  choice  of  the  time  origin,  so  that  for  *VT, 
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2 2 
a 1 V 

q * 


2 2 2 

a - V* 

qy  q y 


a - a2  V v 
aqxy  q x y 


where 

a2  - (AT  )2/12  (24) 

q q 

end  Vx  end  Vy  ere  the  components  of  the  true  velocity. 

2.4  STATISTICAL  PERFORMANCE  OF  THE <*-  P TRACKING  FILTER. 

The  etetieticel  properties  of  the  lest  three  terms  in  (20)  depend  on  the 
noise  reduction  performence  of  the  tracking  filter.  The  trecking  filter 
(in  both  dimensions)  cen  be  considered  es  e single-input  (X,,,),  multiple  out- 
put (XB,  Vfl  Xp)  filter  which  cen  be  shown  to  heve  unbiesed  outputs 
if  Xg,  is  unbiesed,  es  it  is  in  this  cese  (this  is  elso  true  in  the  esse 
where  time-correction  is  used).  Therefore,  the  stetisticel  characteristics 
of  a constant  coefficient,  isotropic «-  8 trecking  filter  are  completely 
expressed  in  terms  of  the  normalised  variance  reduction  ratios, 


where  Ks,  Kv,  and  Kp  are  the  reduction  ratios  for  Xs,  V8,  and 
Xp,  respectively,  and  <?£  is  the  variance  of  the  noise  at  the  input  of 
tne  filter.  Expressed  in  terms  of  the  smoothing  constants,  the  variance 
reduction  ratios  are  [11,  23], 

K - 2 «2  ♦ B(2-3<0  (2* 

a a(4-2a-g) 


!_  2 B 

2 a(4-2'»-p) 

T 


where  T ie  Che  constant  smoothing  interval 


The  variance  reduction  ratios,  however,  do  not  include  the  effects  of  tisw- 
correction  which  is  sisiply  a feedback  loop,  as  shown  in  Figure  1,  in  which 
there  is  an  additional  noise  source  to  account  for  tisw  quantisation  in  the 
tiae-correction  process.  Since  the  errors  are  additive  and  statistically 
independent 


so  that  the  variance  at  the  input  to  the  filter  in  which  time-correction  is 
used  is  o?. 


TRACKING 

FILTER 


FIGURE  1.  FEEDBACK  LOOP  ILLUSTRATING  TIME  CORRECTION 
AND  TIME  QUANTIZATION  ERRORS 


Since  the  introduction  of  time-correction  would  be  expected  to  modify  the 
noise  rejection  performance  of  the  tracking  filter,  the  variance  reduction 
ratios  for  the  case  where  time-correction  is  used  were  calculated  using  the 
standard  z-transform  method  [24],  and  for  this  case 


3ag  ♦ 2B  ♦ fl  AT/T)/a 


2(3/T)  /A 


Kb  - (2a  ♦ aB  ♦ 2B  ♦ B AT/T)/a 


where 


<»  (4-2<»-0)  - 6 (4-4<»-0)AT/T  - 2(e  AT/T)^ 


which  reduce  to  the  reeuite  given  previously  when  AT*0.  The  variance  reduction 
ratios  can  alao  be  derived  using  a matrix  approach  [21,  22]. 

Since  the  three  filter  outputs  are  obtained  from  a common  input,  X,,,,  it 
would  be  expected  that  a nonsero  correlation  would  exist  between  the  various 
outputs  «o  an  additional  item  of  interest,  which  will  be  required  to  deter- 
mine the  statistical  characteristics  of  At,  is  the  covariance  between  the 
velocity  and  the  smoothed  position  which  is  given  by  the  normalised  covari- 
ance ratio, 

1L..  • 0(2c»-0)/(TA).  (36) 


In  the  derivation  of  the  above  equations,  it  was  assumed  that  AT  was  a known 
constant,  yet  in  practice  this  is  not  true.  The  value  actually  used  in  the 
tracking  computations,  ATC,  contains  an  error  resulting  from  time  quanti- 
sation, yet  the  true  value,  ATy,  is,  in  fact,  a random  variable  which  varies 
from  sensor  to  sensor  and  with  the  trajectory  of  the  target  (even  a constant 
rotation  rate  sensor  will  not  produce  a constant  time  interval  unless  the 
target  of  interest  is  stationary).  As  a result,  there  is  some  question  as  to 
what  value  should  actually  be  used  for  AT.  If  the  reference  time  used  by 
the  tracking  filter  for  smoothing  and  prediction  is  taken  as  the  center  of 
the  filter  cycle,  then  this  will  minimise  the  maximum  time-correction  required 
so  that  AT  will  be  restricted  to  the  range  -T/2  to  T/2.  In  order  to  simplify 
the  analysis,  it  will  now  be  assumed  that  AT  is  a constant  which  will  be  set 
equal  to  its  worst-case  value  of  T/2,  since  this  value  results  in  larger  vari- 
ance reduction  ratios  (i.e.,  a higher  noise  variance  at  the  output  of  the 
filter).  As  a result  of  this  assumption  the  statistical  properties  of  all 
terms  in  (20)  are  now  known  with  the  exception  of  <V^. 

The  probability  density  function  for  the  product  term,  <TTe,  can  be  shown  to 
be  an  exponential  integral,  but  since  these  random  variables  are  independent 
with  mean  aero,  the  variance  of  the  product  is  simply  the  product  of  the 
variances.  In  this  case,  the  variances  of  the  last  three  terms  in  (20)  can 
be  obtained  from  the  variance  reduction  ratios  applied  to  o2  so  that,  noting 
that  Va  and  "r^  are  not  independent,  x 


{a2  ♦ V2  a2)  (1  ♦ (AT2  * a2)  K ♦ K 

x x q q v P 

-2AT  (Kvg  ♦ TKV) ) 


where  <;2r  is  the  variance  of  the  x-component  of  Ar  and  similarly  fora2^ 

and  ox yr.  It  ia  implicitly  assumed  in  (37)  that  the  scan-to-scan  differences 
in  the  statistical  properties  of  o2  (and  a2,  ®xy  also)  are  negligible  which 

X y J 


is  a reasonable  assumption  except  at  points  close  to  the  radar  where  the  range 
to  the  target  may  be  on  the  same  order  of  magnitude  as  the  distance  moved  by 
the  target  in  the  smoothing  interval  T.  The  probability  for  any  particular 
search  area  of  interest,  including  the  effects  of  tracking  and  time-correction, 
can,  therefore,  be  obtained  by  integrating  the  bivariate  Gaussian  distribution, 
defined  by  <^r»  tfyr»  and  °xyr»  over  the  region  of  interest. 

The  value  of  the  integral  just  defined  will  be  the  probability  of  a suc- 
cessful fine  correlation  when  the  search  area  for  £r  is  the  same  as  the  area 
over  which  the  integration  is  performed.  The  use  of  the  Gaussian  distri- 
bution requires  the  implicit  assumption  that  the  sensor  measurement  errors 
are  significantly  larger  than  the  errors  introduced  by  finite  precision  compu- 
tations and  time  and  signal  quantization  which  will  be  the  case  in  a well- 
designed  system. 


3. 


NUMERICAL  RESULTS  FOR  CIRCULAR  SEARCH  AREAS 


For  the  purposes  of  this  study  circular  search  areas  will  be  used  even  though 
the  track  datum  deviation  density  assumes  a highly  elliptical  shape  as  the 
range  to  the  target  increases  (resulting  from  the  predominance  of  azimuthal 
errors  at  large  distances).  The  circular  search  area  has  a significant 
practical  advantage  in  that  the  correlation  decision  can  be  based  solely  on 
the  magnitude  of  the  track  datum  deviation,  |&r|,  or  the  square  of  the 
magnitude  in  order  to  avoid  a square-root  computation,  and  the  angular  rela- 
tionship between  the  datum  and  the  track  can  be  ignored.  Thus,  the  circular 
search  area  has  a considerable  computational  advantage  over  noncircular 
search  areas  which  would  impose  a much  higher  computational  burden.  The 
design  of  the  search  area  is  now  reduced  to  the  problem  of  determining  the 
maximum  value  of  the  track  datum  deviation  for  which  correlation  will  be 
allowed.  Since  it  is  desirable  to  maintain  the  same  level  of  tracking  perfor- 
mance throughout  the  entire  coverage  region,  the  size  of  the  search  area  must 
be  a function  of  the  range  from  the  sensor  to  compensate  for  the  increasing 
significance  of  azimuthal  errors.  Acceptable  tracking  performance  requires  a 
high  degree  of  correlation  so  that  the  search  area  must  be  large  enough  to 
include  a significant  portion  of  the  region  in  which  the  desired  data  point 
might  be  found.  To  design  search  areas  to  have  a constant  probability  of 
intercept,  «av  95  percent,  it  is  necessary  to  integrate  the  bivariate  Gaussian 
distribution  over  a circular  region  and  to  present  these  results,  in  terms  of 
the  radius  of  the  circle,  as  a function  of  the  distance  to  the  sensor,  which 
will  achieve  the  desired  level  of  correlation.  To  simplify  the  numerical 
computations  required,  it  would  be  desirable  to  reduce  the  two-dimensional 
integral  to  a single  dimensional  integral,  and  three  means  have'been  found  by 
which  this  can  be  accomplished. 

3.1  EVALUATION  BY  DIRECT  INTEGRATION  OF  f(r). 

The  magnitude,  r,  of  the  radius  vector  for  a correlated  bivariate  Gaussian 
distribution  with  zero  mea.  is  given  by  the  Nakagami  q-distribut ion  [25,26]: 
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(38) 


f<r)  **“  l-f  (t-i)) 

where 

2 2 / 2 2 2 2 
* -aJ  ♦ aJ  ♦ VC*.-®)  ♦ (2co  a ) 

7 * * 7 


(39) 


b 


♦ (2c« 

x 


a 

y 


) 


c “ °xy! 


(40) 

(41) 


end  I0(*)  is  Che  modified  Beesel  function  of  the  first  kind  with  order  sero 
(which  is  essily  celculeted;  e.g. , (27,  28 )X  The  required  probsbility,  using 
(37)  to  cslculste  the  vsrisnces  for  (39)-(41),  is 

P(R)  - f(r)  dr  (42) 


which  csn  be  obtsined  using  one  of  the  msny  numericsl  techniques  for  evsl- 
usting  single  disiensionsl  integrsls  (e.g.,  Csussisn  qusdrsture  [27,  29]  wes 
used  in  this  study). 

3.2  EVALUATION  IN  POLAR  COORDINATES . 

Since  the  rsnge  snd  ssimuth  errors  sre  assumed  to  be  Independent,  equations 
(8)  snd  (9)  csn  be  spplied  using  (37),  sssuming  thst  the  time-qusntisstion 
errors  sre  insignificsnt , to  revert  to  s polsr  coordinste  system  in  which  the 
vsrisbles  of  the  joint  density  function  sre  sepsrsble.  The  limits  on  the 
rsnge  integrstion  ss  s function  of  ssimuth  csn  be  determined  from  the  Lsw 
of  Cosines  (see  Figure  2): 

rmsx  " P cos  0 + /r2  ” ( P sin  6 )2  (43) 

rmin  * P cos  0 “ /r2  - ( P sin  $ )2  (44) 

then,  expressing  the  integrsl  over  rsnge  in  terms  of  the  error  function  l 26 ] 
snd  noting  the  symmetry  in  ssimuth, 
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(45) 


•in  _1(R/p) 

*(*)  • ^•rf((r1M|X-p)/  /Tap)  - «rf ((rBin-p)//2a^)  ) 

•xp<-02/2^) 


which  can  also  be  evaluated  ueing  Gauaaian  quadrature  and  where  erf(*)  it  the 
•tandard  error  function  integral. 


FIGURE  2.  DERIVATION  OF  RANGE  LIMITS  FOR  PROBABILITY  CALCULATIONS 


3.3  EVALUATION  IN  CARTESIAN  COORDINATES. 


1 


The  probability  density  function  of  the  bivariate  Gaussian  distribution, 
f(x,y)  contains  terms  with  the  product  xy  and  so  cannot  be  factored,  but  if 
the  joint  density  is  written  as 

f(x,  y)  ■ f(x)  f(y/x)  (46) 

then  the  integration  can  be  more  readily  accomplished,  because  the  conditional 
density  function  f(y/x)  is  Gaussian  with  moments  (30): 

E(y/x)  ■ c °JL  x 
°x 

E(y^/x)  ■ (l-c^)  ♦ c2  x2 

•i 

where  c is  the  correlation  coefficient  and,  using  this,  the  variables  in  (46) 
are  separable.  Performing  the  integration  over  y first  yields  an  integral 
of  the  same  form  as  (43)  in  which  the  result  of  the  first  integration  is  the 
difference  of  two  error  functions  in  which  the  limits  of  integration  on  y 
are  expressed  in  terms  of  x and  which  can  also  be  evaluated  using  Gaussian 
quadrature. 

Since  a considerable  amount  of  computation  was  required  in  this  study,  the 
three  methods  discussed  above  were  developed  to  provide  computationally  effi- 
cient techniques  for  evaluation  of  the  integrals  required.  Other,  more 
generalised,  techniques  have  also  been  developed  for  integration  of  the 
multivariate  Gaussian  distribution  of  higher  dimension  [31) . 

3.4  DISCUSSION  OF  NUMERICAL  RESULTS. 

For  computational  purposes  the  following  parameter  values  have  been  used 
(corresponding  to  those  for  the  enroute  air  traffic  control  system): 

°pm  0.125  nmi,  o £■  0.263°,  Vx  ■ Vy  ■ 600  knots,  AT„  “0.5  a.ATj  ■ 3.0  s, 

T ■ 6.0  s,  a ■ 0.3125,  and  0 • 0.046875.  The  values  chosen  for  Vx,  Vy,  and 
ATT  were  selected  es  the  worst  case  values  which  would  be  observed  in  normal 
situations.  Using  these  parameter  values,  the  contribution  of  timing  errors 
in  (37)  is  negligible,  with  the  result  that  the  variance  is  about  28.1  percent 
larger  than  the  variance  of  the  measuresmnt  error  alone,  and  this  is  due  mainly 
to  the  contribution  of  the  predicted  position  Kp.  To  obtain  the  radius  of 
the  circle  required  to  assure  a specified  level  of  correlation,  in  this  case 
95  percent,  a numerical  interpolation  routine  was  used  to  find  the  radius 
as  a function  of  the  desired  probability.  The  numerical  interpolation  was 
required  because  it  is  considerably  easier  to  find  the  probability  as  a 
function  of  the  radius,  rather  than  vice  versa.  The  results  of  these  compu- 
tations are  given  by  the  lower  line,  corresponding  to  a xero  offset,  in 
Figure  3.  Note  that  the  results  are  independent  of  axismth  and  depend  only 
on  the  distance  to  the  sensor.  To  insure  computational  accuracy,  both 
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ZUO  BIAS  +1.5 


FIGURE  3.  SEARCH  AREA  RADIUS  FOR  95-PERCEWT  CORRELATION 


msthods  1 and  2 above  vara  uaad  Co  provide  a check  on  Che  reaulca.  In 
Che  high-probability  region*,  which  are  Che  only  regions  of  inCeresC  in  Chi* 
case,  Che  reaulca  of  Che  Cvo  vet hod*  cypically  differed  by  abouC  2 x 10“^ , 
which  is  More  chan  sufficienc  accuracy  for  Che  purposes  of  chis  sCudy.  A* 
would  be  expecCed,  Che  radius  of  Che  search  area  is  relacively  constant  at 
short  distances  and  becomes  linearly  proportional  to  distance  at  larger  ranges 
in  Che  region  where  asimuthal  errors  predominate. 

For  convenience,  the  results  for  Che  case  where  Che  mean  deviations  are  not 
saro  are  also  given  in  Figure  3.  Nonxero  expecCed  values  arise  when  Che 
predicted  position  is  offset  from  Che  true  target  position  which  usually 
results  from  a maneuvering  target  so  that  the  aasuswd  target  dynamics  are 
not  correct  and  the  first  three  terms  of  (19)  will  be  nonxero.  The  results 
in  this  case  were  found  to  depend  on  the  sMgnitude  of  the  offset,  do,  and 
on  the  angle  between  the  line  from  the  target  to  the  origin  and  the  line  from 
the  target  to  the  point  of  offset,  as  illustrated  by  the  angle0in  Figure  4. 

The  results  were  also  found  to  be  periodic  in^with  the  maximum  and  minimum 
values  occurring  at  0*  and  90*  so  that  only  these  two  cases  need  be  considered, 
since  only  the  radius  of  the  search  area  is  of  importance.  The  fact  that  the 
auixiaxim  and  minimum  of  the  search  area  radius  occur  at  either  0*  or  90*  is  due 
to  the  predominance  of  either  range  or  asimuthal  errors,  with  the  crossover 
point  at  the  range  where  the  errors  are  equal  in  magnitude.  Using  (3)  to  (7), 
the  x and  y variances  will  be  equal  and  the  covariance  xero  if 


(49) 


or 

P - op/"*, 


(30) 


and  using  the  standard  deviations  given  previously,  the  crossover  point  ia 
27.2  nmi  as  illustrated  in  Figure  3.  The  numerical  computations  for  the 
biased  case  were  performed  by  modifying  the  second  and  third  methods  dis- 
cussed previously  to  account  for  nonxero  means  in  the  x and  y track  datum 
deviations  with  the  offset,  do,  as  the  magnitude  of  the  means: 


■(Ax)  ♦ E2(Ay). 


(51) 


For  computational  purposes,  the  first  method  used  previously  is  impractical 
in  this  case  because  the  probability  density  function  is  given  by  an  infinite 
series  in  which  the  individual  terms  are  expressed  as  modified  Bessel  function* 
125,  32). 


Since  it  is  the  objective  in  the  design  of  a search  area  to  achieve  a high 
probability  of  correlation,  it  is  apparent  that  the  search  area  must  cover  a 
significant  portion  of  the  probability  density  function  to  achieve  this  goal. 
It  would  be  expected,  therefore,  that  the  radius  of  the  search  area  in  the 
biaaad  case  would  have  to  be  sufficiently  large  so  as  to  compensate  for  most 
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of  the  bias,  and  tha  raaulta  in  Figure  3 show  thia  Co  ba  cha  case;  however, 
it  ia  not  nacaaaary  to  compensate  for  Cha  anCira  amount  of  Cha  bias  as  shown 
by  Cha  rasulcs,  alao  in  Figura  3,  in  which  1.5  nai  was  aiaply  added  Co  cha 
radiua  required  in  Che  unbiaaed  caae.  When  a biaa  occurs,  ic  ia  likely  chac 
Che  Cargec  ia  maneuvering,  in  which  case  differenc  smooching  constants  aay  be 
used.  Also  shown  in  Figure  3 are  Che  resulca  for  che  caae  where  Che  biaa  was 
1.5  nai  and  che  a and  3 smooching  conscancs  were  0.5  and  0.15625,  raspecCively, 
and  chase  resulca  agree  fairly  closely  wich  Chose  for  Che  case  where  1.5  nmi 
was  added  co  Che  radius  requirad  in  Che  unbiaaed  caaa  (which  were  obcained 
wich  differenc  ssnoching  conscancs).  For  praccical  purposas,  Cherefore,  Che 
radius  of  a circular  search  area  required  for  a constant  probability  of  cor- 
relation can  be  taken  as  che  sum  of  Che  radius  required  in  Che  unbiased  case 
plus  che  magnitude  of  any  biaa  which  may  be  present,  and  this  ia  especially 
Crua  in  che  case  where  larger  sanothing  constants  are  used  when  a biaa  is 
detact ed. 


FIGURE  4.  ILLUSTRATION  OF  OFFSET  SEARCH  AREA  FOR  NONZERO  BIAS 
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3.5  TRACKING  BIASES  FOR  MANEUVERING  TARGETS. 


It  vaa  seated  previously  that  one  of  the  functions  of  the  correlation  task  is 
to  choose  between  alternative  sets  of  smoothing  constants,  and  this  function 
can  be  achieved  by  dividing  the  correlation  region  into  subareas,  each  of 
which  is  associated  with  a particular  set  of  smoothing  constants.  For  the 
purposes  of  this  study,  the  correlation  region  will  be  divided  into  large 
and  small  search  areas  as  illustrated  in  Figure  5.  The  sixe  of  the  ssmII 
search  area,  in  which  smoothing  constants  applicable  to  straight-line  tracks 
are  used,  is  determined  by  the  results  in  Figure  3 for  the  unbiased  case, 
while  the  sixe  of  the  large  search  area,  in  which  smoothing  constants  appli- 
cable to  turning  tracks  are  used,  will  be  taken  as  the  sum  of  the  radius 
of  the  small  search  area  plus  the  magnitude  of  any  bias  which  may  be  present. 
Therefore,  in  order  to  determine  the  sixe  of  the  large  search  area  it  is 
necessary  to  determine  the  magnitude  of  the  bias  which  will  result  from  real- 
istic maneuvers.  By  determining  the  maximum  bias  which  can  be  realistically 
be  expected  to  develop,  it  will  be  possible  to  design  the  large  search  area 
on  a worst  case  basis  so  that  correlation  will  be  possible  for  all  realistic 
maneuvers. 


Implicit  in  the  equation  for  the  predicted  position  is  the  assumption  of  a 
straight-line  constant  velocity  trajectory,  and  any  deviation  from  this  ideal- 
ixed  trajectory  is  likely  to  cause  a bias  error.  Since  the  significance  of 
the  bias  will  vary  in  different  situations,  the  comparison  of  the  bias  will 
be  in  terms  of  the  relative  bias  as  normalised  with  respect  to  the  prediction 
distance,  VT.  For  example,  if  the  target  accelerates  in  a straight-line 
trajectory,  then  the  relative  or  normalised  bias  is 


VjaT2/VT  - SaT/V. 


(52) 


A PREDICTED  POSITION 
J SMALL  SEARCH  AREA 
^ LARGE  SEARCH  AREA 
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ILLUSTRATION  OF  LARGE  AND  SMALL  SEARCH  AREAS 
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Assuming  worst  esse  values  of  100  knots/minute  and  100  knots  for  acceleration 
and  velocity,  respectively,  for  a 6-second  sampling  time,  the  normalized  bias 
is  0.05,  but  in  most  cases  it  will  be  less  than  this  and  negligible  for  prac- 
tical purposes.  A climbing  or  descending  target  can  also  introduce  a bias 
in  the  predicted  position.  If  a target  is  moving  at  a constant  velocity  and 
then  begins  a vertical  transition  at  a constant  vertical  velocity,  Vv,  while 
maintaining  the  same  airspeed,  then  the  true  ground  speed  will  decrease, 
resulting  in  a normalized  bias  of 

dn  - 1 - 1l  - (Vv/V)2  (53) 

and  choosing  worst  case  values  of  4,000  feet/minute  and  100  knots  for  Vv  and  V, 
respectively,  the  normalized  bias  in  this  case  is  0.085,  which  is  also  negligi- 
ble for  practical  purposes. 

Both  of  the  biases  just  discussed  are  relatively  small  and  occur  in  the  line 
of  the  trajectory.  For  a turning  trajectory,  however,  substantially  larger 
biases  can  develop  which  are  not  along  the  straight-line  projection  of  the 
velocity  vector.  For  a target  moving  at  a constant  speed  in  a circular 
path  in  which  the  heading  changes  at  a constant  rate.co,  it  can  be  shown,  by 
integration  of  the  parametric  equations  of  motion,  that  the  normalized  bias 
at  time  T is 

d,,  - Vl  ♦ (sin(‘^T/2)/(<uT/2))2  - fcJTsinWT)  (54) 

which  is  plotted  in  Figure  6.  The  results  show  that  the  normalized  bias 
increases  with  the  total  heading  change, tuT,  to  a maximum  of  1.26  and  then 
asymptotically  approaches  one  in  an  oscillatory  manner.  The  reason  the  normal- 
ized bias  approaches  one  is  illustrated  in  Figure  7.  As  the  heading  change 
increases,  with  the  target  repeatedly  traveling  a circular  trajectory,  a point 
is  soon  reached  at  which  the  diameter  of  the  trajectory  is  much  less  than  the 
linear  projection,  hence  the  asymtotic  value  of  one.  Since  a typical  turning 
rate  is  on  the  order  of  3*  per  second  or  less,  it  is  apparent  that  a signifi- 
cant bias  would  not  develop  until  several  scans  had  elapsed  during  which  no 
data  were  received  for  the  track.  If  a target  maneuver  takes  several  scans, 
as  it  usually  does,  then  the  effect  of  tracking  on  the  bias  must  also  be 
cons idered . 

It  has  been  shown  chat  of  the  maneuvers  which  have  been  considered,  the  bias 
observed  during  a turn  is  the  largest,  however,  this  is  without  considering 
the  effect  of  tracking.  For  simplicity  in  analyzing  the  effect  of  tracking  it 
would  be  desirable  to  restrict  consideration  to  turning  maneuvers.  The  justi- 
fication for  this  simplifying  restriction  is  that  for  nonturning  maneuvers, 
only  the  constant  velocity  assumption  used  in  the  tracking  prediction  is 
violated,  while  a turn  represents  a violation  of  both  the  straight-line  and 
constant  velocity  assumptions.  In  addition,  a turn  represents  a time-varying 
acceleration  (in  Cartesian  coordinates),  while  the  other  maneuvers  considered 
represent  a constant  acceleration.  It  is  concluded,  therefore,  that  if  the 
bias  which  is  used  to  determine  the  size  of  the  large  search  area  is  based 
on  a turn,  which  represents  a worst  case  situation,  then  if  it  is  possible 
to  track  turning  targets  with  this  large  search  area  size,  it  will  also  be 
possible  to  track  other  maneuvers  which  produce  a smaller  bias. 
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PREDICTED  POSITION 


FIGURE  7.  ILLUSTRATION  OF  BIAS  FOR  CIRCULAR  TRAJECTORIES 


Unless  the  tracking  function  is  exactly  able  to  compensate  for  the  bias  as 
it  occurs,  which  is  not  possible,  then  it  would  be  expected  that  the  effect  of 
the  bias  on  each  scan  would  be  cumulative  and  would  increase  until  a steady - 
state  value  had  been  reached.  The  magnitude  of  the  steady-state  bias  in 
position  has  been  determined  by  calculating  the  magnitude  of  the  sinusoidal 
response  of  the  tracking  filter  transfer  function,  which  would  correspond  to  a 
circular  trajectory,  and  after  normalization 


d 

n 


4 sin2(ct>T/2)/a/r 

{[(a+6-2-SAT/T)  + (2-a+BAT/T)  coscJT]2 

2 

+ [(a-BAT/T)  sinaJT]  }** 


(55) 


where  a and  8 are  the  smoothing  constants  of  the  tracking  filter  and  cjT  is 
the  heading  change  per  scan  (rather  than  the  total  heading  change  as  pre- 
viously used).  The  normalized  bias  calculated,  using  (55),  agrees  with  the 
results  obtained  in  another  study  [9]  and  is  given  in  Figure  8 for  the  two 
parameter  sets  used  previously  and  for  AT  * t T/2  and  0.  The  results  in 
Figure  8 show  that  the  peak  bias  errors  developed  during  tracking  of  a 
maneuver,  which  occur  in  the  coT-region  of  practical  interest,  can  be  two  to 
three  times  larger  than  the  projected  distance  for  a straight-line  trajectory, 
VT,  and  also  that  the  process  of  time-correction  has  little  effect  in  this 
case.  Clearly,  the  bias  developed  during  the  tracking  of  a sustained  maneuver 
is  the  largest  bias  which  must  be  accommodated  in  the  design  of  the  large 
search  area.  Presumably,  by  the  time  a bias  of  this  magnitude  had  developed, 
the  tracker  would  have  switched  to  the  large  search  area  smoothing  constants 
so  that  the  maximum  size  of  the  large  search  area  for  purposes  of  correlation 
could  be  taken  as  2. A times  the  velocity-time  product  plus  the  small  search 
area  radius.  Since  the  velocity  which  must  be  used  for  this  purpose,  how- 
ever, is  the  estimated  velocity,  given  by  (11),  it  will  be  necessary  to  use 
a different  multiplier  than  that  indicated  for  the  position  bias  in  Figure  8, 
because  the  estimated  velocity  during  a maneuver  becomes  significantly  less 
than  the  true  velocity. 

The  magnitude  of  the  steady-state  bias  in  the  velocity  can  be  determined  in 
the  same  manner  as  that  used  for  the  position;  i.e.,  the  magnitude  of  the 
sinusoidal  response  of  the  tracking  filter  transfer  function  for  velocity  will 
be  the  steady-state  speed.  Using  this  approach,  the  ratio  of  the  estimated 
speed  to  the  true  speed  is 

{[  (3-a-8+BAT/T)  + (2a+6-A-2&AT/T)cos<jr+(l-a+eAT/T)cos2<jr]2lJi 

+[  ( 2a*  3-2-2  8 AT/T ) s in£dT+  ( l-a-8  AT/T)  s in2cJT  ]2  ! 

■ -a-  , p_ __  (56) 

lVTl  ^ {C(a+e-2-eAT/T)  + (2-a+8AT/T)cosafl']2  + [(<*"BAT/T)sin  T]2  } ** 

which  is  also  plotted  in  Figure  8.  As  the  results  show,  the  bias  in  velocity 
can  be  significant,  especially  for  the  small  search  area  smoothing  constants 
given  previously.  To  check  the  accuracy  of  this  equation,  a statistical  simu- 
lation of  the  tracking  algorithm  was  performed  and  the  simulation  results,  in 
terms  of  the  ratio  of  the  estimated  speed  to  the  true  speed,  agreed  very  well 
with  the  ratio  obtained  using  (56). 

If  the  estimated  velocity,  which  in  practical  situations  is  the  only  velocity 
information  available,  is  to  be  used  for  determination  of  the  large  search 
area  size,  then  in  order  to  compensate  for  the  velocity  bias,  the  estimated 
velocity  must  be  divided  by  the  normalized  velocity  bias  factor  in  order  to 
obtain  an  estimate  of  the  true  velocity.  Therefore,  if  the  product  of  the 
filter  period  and  the  estimated  velocity  is  to  be  used  to  calculate  the 
maximum  position  bias  in  order  to  specify  the  size  of  the  large  search  area, 
then  this  product  must  be  multiplied  by  the  ratio  of  the  normalized  position 
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FIGURE  8.  NORMALIZED  POSITION  AND  VELOCITY  BIASES 
DURING  MANEUVERS  WITH  TRACKING 
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bias  to  the  normal iaed  velocity  biaa  in  order  to  compensate  for  both  the 
position  and  velocity  biases.  The  bias  ratios  for  the  two  sets  of  smoothing 
parameters  considered  previously  are  also  given  in  Figure  8,  and  it  can  be 
seen  that  there  ia  an  almost  linear  relationship  on  a logarithmic  scale  and 
alao  that  the  effect  of  time-correction  is  eliminated  in  these  ratios.  Since 
one  of  the  main  purposes  of  t isw-correct ion  is  to  allow  the  use  of  sensors 
which  operate  with  a different  period  than  that  of  the  tracking  filter,  then 
the  correlation  process  must  be  designed  to  accommodate  such  a situation. 

Conaidering  that  a typical  turn  rate  is  usually  less  than  3*  per  second  and 

that  the  maximum  period  of  rotation  for  ATCRBS  sensors  is  12  seconds,  then 

the  maximum  biaa  ratio  required  for  the  large  search  area  is  four  (from 

Figure  8 with  a>T-36),  which  ia  considerably  larger  than  the  maximum  value  of 
2.4  found  previously  by  considering  only  the  position  bias.  Therefore,  in 
order  to  guarantee  correlation  in  the  large  search  area  with  at  least  the  same 
probability  aa  for  the  small  search  area,  the  sise  of  the  large  search  area 
should  be  taken  as  the  sum  of  the  small  search  area  sice  plus  four  times  the 
magnitude  of  the  product  of  the  scan  time  and  the  estimated  speed. 


4. 


CONCLUSIONS 


The  results  just  presented  provide  the  basis  for  the  design  of  correlation 
regions  on  a conatant-probability-of-intercept  basis.  In  many  cases,  two  cor- 
relation regions  are  required:  a small  search  area  used  for  straight-line  tracks 
and  a large  search  area  used  for  maneuvering  targets  with  different  smoothing 
constants  for  each  search  area.  The  sice  of  the  small  search  area  is  deter- 
mined primarily  by  the  measurement  errors  in  the  sensor  and  the  tracking 
performance  for  straight-line  tracks.  At  distances  close  to  the  sensor  a 
constant  radius  is  used  for  the  small  search  area,  while  as  the  distance  to 
the  sensor  increases  the  search  area  radius  becomes  linearly  proportional  to 
distance.  The  radius  of  the  large  search  area  is  equal  to  the  sum  of  the 
radius  of  the  small  search  area  plus  the  magnitude  of  the  worst  case  bias 
which  could  reasonably  be  expected  to  develop  in  a sustained  maneuver.  The 
specification  of  the  sise  of  the  correlation  region  in  this  manner  assures  a 
constant  level  of  performance  with  regard  to  the  probability  of  correct 
correlation  and,  in  addition,  results  in  better  discrimination  against  false 
targets  and  improved  performance  in  turn  detection  for  purposes  of  tracking. 

The  use  of  a dynamically  varying  correlation  region,  as  discussed  in  this 
study,  was,  in  fact,  found  to  yield  improved  tracking  performance  as  compared 
to  the  performance  obtained  using  fixed  radius  correlation  regions  |33). 

Although  the  numerical  results  of  the  present  study  apply  to  a specific  a-t* 
tracking  filter,  the  technique  developed  could  easily  be  applied  to  other 
tracking  filters. 
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